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I. INTRODUCTION 



SN1987A neutrino events [l| prompted many dedicated 
analyses. Even if the number of detected neutrinos is lim- 
ited, these analyses provide interesting limits on neutrino 
properties and clues on the core collapse mechanism. 



The question of which likelihood should be adopted 
for supernova neutrino data analysis will become crucial 
after the next galactic supernova, when a much larger 
number of neutrino events will be collected. These events 
will carry information on neutrino properties mixed with 
information about supernova properties, so that we will 
need to jointly study their energy, time and direction 
distributions to try to extract all the relevant pieces of 
information. Therefore it will remain unpractical to bin 
the events and an event-by-event likelihood will remain 
the best tool for data analysis. 



We present here the likelihood that should be ideally 
adopted for supernova neutrino data analysis. Our like- 
lihood is more general than those already present in the 
literature @, Moreover, we resolve discrepancies 

in the previous literature, numerically relevant already in 
the concrete case of SN1987A data. We argue, in particu- 
lar, that the analysis of SN1987A neutrino data by Lamb 
and Loredo || (LL), quoted since 2004 in the summary 
table of the Particle Data Group @, uses a likelihood 
that incorrectly biases the analysis in favor of low energy 
events. We here present the correct likelihood, generaliz- 
ing the 'traditional' form, advocated, e.g., by Jegerlehner, 
Neubig and Raffelt 



The structure of this paper is the following. In Sect. ITT1 
we derive the general form of the likelihood. The applica- 
tion to a specific case of interest is discussed in Sect. IIIII 
Finally, in Sect. lIVl we compare our likelihood with other 
forms adopted for the analysis of SN1987A neutrinos, 
showing how the fitted parameters got biased. 



II. DERIVATION OF THE LIKELIHOOD 



A. General form of the likelihood 



We write the expected event number in the i-th bin as: 



= dti ebb 



dN 
dt <ix 



{pi > x i ) j 



(1) 



where ti represents the time coordinate, while x^ indi- 
cates the set of all other observables (energy, position, 
direction, etc.) which define the properties of the i-bin. 
We suppose that the bin sizes dti dx.i are infinitesimally 
small so that the condition rii <C 1 holds true: therefore 
the probability that multiple events are collected in one 
bin is negligible and, thus, observing N cv events corre- 
sponds to N cv bins with 1 event, and all other bins with 
events. 

According to Poissonian statistics (see e.g., Appendix 
A of 0) the associated likelihood is: 



C = cxp 



E 



n 



(2) 



where the sum in the exponent runs over all A^m bins 
and gives the total number of expected events, while the 
product runs over all iV ev observed events. As usual, one 
can convert this into a \ 2 distribution as C — e~ x I 2 . 



B. Distinguishing between signal and background 

Let us consider the case when the detected events are 
due to a signal S, reprocessed in the detector through a 
response function 1Z, and to a known (measured) back- 
ground process B. We have: 



dN 
dtdx 



(i,x)= B(t,x) + J dt'dx'S{t',x')n(t',x',t,x) 



(3) 



The second term in the r.h.s. takes into account that a 
signal produced at the time t' and with coordinates x', 
due to detector response, could be observed with a proba- 
bility lZ(t' , x', t, x) at a different time t and coordinate x. 
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By integrating over all possible detection times and co- 
ordinates, we introduce the general form of the detection 
efficiency: 



T](t',x') = J dt dxK(t\x',t,x) 



(4) 



The efficiency obeys the condition < r) < 1, if we de- 
scribe a situation when the events can be lost. By fac- 
toring out rj we define the smearing (or error) function Q 

Q(t', x', t, x) = H(t', x', t, x)/r)(t\ x') (5) 

normalized to unity: 



J dt dx g(t',x',t,x) = 1. 



(6) 



The background B, the efficiency r\ and the smearing 
Q describe the experimental apparatus. Assuming that 
they are known, we can use an experimental result to 
learn on the signal S, by the study of the likelihood func- 
tion of Eq. © together with ©. 



C. Simplifications 

In the case of interest, it is possible to further simplify 
the problem by relying on the following assumptions: 

(i) We assume that the response function factorizes in 
the time and in the coordinates as follows 



Tl(t', x', t, x) = r(t',t)K(x', x). 



(7) 



We introduce the time-independent efficiency in the ob- 
servables Jj(x'), defined in analogy to Eq. j4|): 



r/(x') = J dx K(x',x), 



(8) 



and the smearing function defined in analogy with 
Eq. ©: 



S(x',x) = ft(x',x)Mx'). 
Again, it is normalized to unity: 



J dx G{x',x) = 1. 



(9) 



(10) 



We will discuss later the specific form of these expressions 
for SN1987A. 

(ii) If the time t is measured with negligible error, we 
have 



r(t',t) = S(t-t'), 



(11) 



possibly multiplied by a window function w(t) to account 
for the dead time r after an event, due to supernova or 
to background (for example, a muon), has been recorded. 
Concerning SN1987A data, only the relative time be- 
tween events of the detectors was measured precisely; one 



needs to take into account the uncertainty in the absolute 
time of the Kamiokande-II and Baksan events. 

(Hi) We can finally assume that the background does 
not depend on the time, namely 



B(t,x) = B(x) 



(12) 



possibly, multiplied by w(t) to take into account for the 
absence of any events, including those due to background, 
during dead time. Eq. (|12|) implies that the background 
can be measured in the period when the signal is absent 
(as for SN1987A). 
With these assumptions, Eq. ([3]) simplifies to: 



dN 
dtdx 



(t,x) = B(x) + J dx'g(x',x)r](x')S(t,x') 



(13) 



Then, assuming that the N av events Xj have been mea- 
sured at time U, the likelihood in Eq. @ becomes: 

£ _ — J dtdxB(x) — J dtdx'r)(x')S(t,x ) 

U^l [B(xi) + Jdx' g(x',x i )ri(x')S(t i ,x')] dUdxi, 

where, in the exponent, we replaced the sum over all in- 
finitesimal bins with an integral and used (|10p . By drop- 
ping constant factors, that are irrelevant for estimating 
the parameters that control the theoretical expression of 
the signal rate S, and replacing x' with x, we get 



C 



- J dt dx 7)(x)S(t,x) 

nti i B (^) + I dx £(*, Xi)T?(x)S(f,, x)] 



(15) 



This form of the likelihood is general enough for the pur- 
pose of analyzing SN1987A neutrinos. Moreover, this is 
a generalization of the likelihood advocated in [8[ for the 
study of radioactive decays, when the time of occurence 
of each event is measured. 

As we already discussed, the dead time can be taken 
into account by extending the time integral in the expo- 
nent only to the time when the detector is on, thereby re- 
moving the time intervals where data taking was stopped 
after each candidate signal event. As long as r is small 
enough, one can equivalently take into account the dead 
time due to background events by multiplying the inte- 
grand in the exponent of (|15[) by the average live-time 
fraction, 1 — rb^, where is the time-averaged back- 
ground event rate. Compare it with the discussion of |{|, 
further elaborated in 



III. APPLICATION TO KAMIOKANDE-II 

In order to specify the general formulae, we choose a 
concrete and important example: we discuss the likeli- 
hood for the water Cerenkov detector Kamiokande-II. 



A. Generalities 



In this subsection, we collect some useful definitions. 
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The variables that characterize an event are: 

x.j = {Ei (energy) , n, L (direction) , r*i (position) }. 

(16) 

For concreteness, we consider events resulting from the 
reaction v e p — > ne + when a positron is detected through 
its Cerenkov light; similar considerations apply to the 
elastic scattering reaction or the charged current reac- 
tions with nuclei. 

In the construction of the likelihood 3 different direc- 
tions are relevant: the direction n* of SN1987A; the re- 
constructed direction hi of each event; the true direction 
h of the positrons produced by the detection process. For 
each event, the first 2 directions are fixed, while we have 
to integrate on the true direction of the positron, taking 
into account the detector response and the reconstructed 
event direction, as described in Eq. (fT5j) . To do this, it 
is convenient to use an "event-centric" system in which: 

1) The reconstructed positron direction is along the z 
axis, 

fk = (0,0,1). (17) 

2) The true positron direction is in the generic direction: 

n= (sin 9 cos p, sin 9 sin ip, cos 9), (18) 

so that 



cos( 



(19) 



Thus, 9 is the opening angle around the reconstructed 
direction and p is the azimuthal angle. The experimental 
collaborations usually quote the error on the angle S9i 
between the true and the reconstructed direction for each 
bin (rather than the the error on the direction versor Shi 
itself). 

3) Finally, we have the versor h* pointing in the direction 
of the supernova. This, without loosing in generality, can 
be chosen in the plane x — z: 



h„ = (sin Si, 0, cos#;), 



so that 



cosf 



Tlfy 71/% • 



(20) 



(21) 



B. Smearing function 

In the simplest approximation, we can describe the 
smearing function by assuming that it factorizes accord- 
ing to: 

S(x,Xi) =G 1 (E-E i ,* 1 )G 2 (h-h i , a 2 )G 3 (f-n,a 3 ) (22) 
where we denote by 

exp(— x 2 /2a 2 ) 



a standard Gaussian in n dimensions. We include a nor- 
malization factor N n to describe the presence of phys- 
ical boundaries, like e.g., the fact that 9 £ [0,7r] and 
ip £ [0, 2tt] when we integrate over the possible directions 
of h. 

The quantities (71,2.3 & re functions of the variables of 
Eq. HH). We identify 0^2,3 in the point E — E^h = 
hi,r = ?i with the error for the i-th event quoted by the 
experimental collaborations, e.g., 

ai(Ei,hi,fi) = SEt 
a 2 {E i ,n i ,f i ) = Srii 
a 3 (Ei,ni,ri) = Sr t (24) 

thus, we further approximate the smearing function: 

£(x,x 4 ) » G 1 (E-E i ,SE i )G 2 (h-h i ,Sn i )G 3 (f f -f f i ,6r i ) 

(25) 

in the vicinity of each point where the likelihood should 
be evaluated — see Eq. (JT3J) - 1 



C. Remarks on the angular distribution 

The expressions for energy and position smearing func- 
tions are essentially standard and do not need particular 
attention. The angular distribution requires, instead, a 
more detailed discussion. 

First, we discuss the connection between the error Srii 
to be inserted into Eq. (I24|) and the 1-sigma error on angle 
59i indicated by the experimental collaborations. In the 
Gaussian assumption, the function G 2 can be written as: 

G 2 in—h; L1 5rii)dh — — ^ x dc—^-(c,8ni) (26) 

where the function dp/dc, given by: 



dp 
dc 



N 2 Sn; 



■ exp 



(27) 



describes the distribution of the angle c = cos 9 — hih 
between the true and the reconstructed direction; the 
azimuthal angle p is uniformly distributed as appropriate 
for an unbiased detector. The normalization factor N 2 
can be explicitely calculated 



N 2 = 1 - cxp(-2/<5^ 2 ), 



(28) 



G n (x, cr) 



A/„(VW)« ' 



(23) 



1 An (arguably) more refined approximation that takes into ac- 
count the Poisson nature of the photoelectron detection can be 
obtained by multiplying the constant errors SEi, Srii and Sri, by 
' Ei/E. We note, in this respect, that the normalization con- 
dition in Eq. (1 1 P is obeyed even in the general case, when the 
functions 01,2,3 vary with the true coordinates of positron, since 
one integrates over the reconstructed coordinates. 
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and it is very close to one for the typical case 5m <C 1. 
We calculate 5n,i by requiring that: 



,w, 



d,9 



dp 
d9 [ 



,6m) = 0.683, 



(29) 



for the 1 sigma error 66 i corresponds to the w 0.683 con- 
fidence level. For small 6m, we get easily: 



6m ^ 0.660 69, x (1 - 69 2 /24). 



(30) 



Typically the first term provides an adequate approxima- 
tion for the quantity we search, 6m- 

Now, we discuss a possible improvement of the Gaus- 
sian assumption for the distribution G^- Experimen- 
tal investigations of the Super-Kamiokande collabora- 
tion [ll[ have shown that the tails of the angular distri- 
bution fall slower than exp(— cte • 9 2 ) and resemble more 
closely exp(— cte • 9); compare also with App. C of \L2\ . 
This suggests to release the Gaussian approximation for 
the distribution on the directions and to replace G2 —> 
exp(— \n— m\/6ni). Thus, the distribution over the cosine 
becomes: 



dp 
dc 



1 



No 6n; 



■ exp 



6m 



(31) 



where the proportionality constant 



2.5 

i 2.0 

1 1.5 
p, 

I 1.0 

§ 

m 0.0 




n/4 
8 [radians! 



71/2 




FIG. 1: Angular distribution (top) and related cumulative 
distribution (bottom) as a function of the angle 6. The con- 
tinuous curves show Eq. (|2T[) and the dashed curves show 
Eq. (|5I). We here assumed 56 = 83° A /MeV/£, as appro- 
priate for Super-Kamiokande [Til ] and a typical energy of an 
elastic scattering event, E — 15 MeV, so that 56 = 0.374. 



JV a = 1 - (1 + 2/Sth) exp(-2/6m), (32) 

is again close to one for small 6m. By imposing the 
condition ([29]) and considering again the limit of small 
6m we calculate the new expression for 6m, obtaining: 

6m ^ 0.424 69i x (1 - 5Q\ /24). (33) 

where, as in Eq. (|30|) . the first term is typically sufficient. 

The two distributions are depicted in Fig. [1] for a spe- 
cific value of 59 i. It is worthwhile to note various features 
of Eqs. ([27]) and (|3i"|) : 

1. When considered as functions of the direction we 
see that they both depend only \h— fii\ and have a 
maximum for h = m, as it should be. 

2. It is easy to treat them analytically, which is a wel- 
come property to use them in a likelihood. 

3. It is simple to study their limit for small 6m by 
replacing sm9 — ► 9 and cos 9 — > 1 — 9 2 /2, which 
makes their analytical treatment even simpler. 2 



2 When only the second replacement is done, E q. H27H coincides 
with the form commonly used in the literature llfi and Eq. (I3U 
practically coincides with the form given in 1131 • In fact, the 
exponential term in Eq. CI of |12| | can be neglected in comparison 
to the linear term x for all relevant energies. 



4. For a fixed 59 i, we see that 5m is smaller in the 
second case; thus, the maximum at cos9 = is 
higher in the second case. 

5. The most probable angle is 9 ~ 5m in both cases; 
thus it is smaller in the second case. 

A choice between Eq. (|2T|) and Eq. (f3"Tj) (or other rea- 
sonable approximations) is not critical for the analysis of 
SN1987A in view of the limited event sample. However, 
the use of an appropriate distribution is potentially im- 
portant for the analysis of the elastic scattering events 
from a future supernova in a water Cerenkov detector. 



D. From the idealized to the actual likelihood 

We are now in the position to provide a concrete and 
useful expression for the likelihood that takes into ac- 
count the reported information on the data and on the 
detector response. 

We recall that: 
1) The signal is expected to be uniformly distributed in- 
side the detector. The angular dependence of the sig- 
nal arises from the angular distribution of positrons pro- 
duced by v e p — > ne + . This can be expressed as a function 
of the angle nn, between the direction h of the produced 
positron and the direction n* of the SN1987A. We can 
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thus replace in Eq. (fT5|) : 

S(t, x )- 



S(t, E, nn«) 
2nV 



(34) 



where V is the volume of the detector, the product nn* 
can be expressed through Eqs.f IT5]) and (|2T))) as: 



cos W.; cos ( 



sin sin a cos 95, 



(35) 



and the factor 2ir accounts for the fact that positron di- 
rections are uniformly distributed with respect to rota- 
tions around n*. 

2) We assume that the background does not depend on 
time, direction and position. We indicate the total back- 
ground counting rate as a function of the energy with 
B{E) and we replace in Eq. (JTE 



B(x,) 



B(E) 

4:TTV ' 



(36) 



3) The average efficiency of the detector as a function of 
the energy fj(E) is known. We assume that the efficiency 
does not depend on time, position and direction, so that 
we can replace in Eq. (|15jl: 

(37) 



rj(E,n,r) n(E) ; 



4) The errors on the energy SEi and on the direction 
59i in the neighbourhood of the given datum are known. 
We additionally indicate by Sri the value of the error on 
the position, on which we have only limited information, 
5r l ~ 1 m at 10 MeV. 

At this point, we have all the elements to write the 
concrete form of the likelihood. Integrating away the 
Gaussian on the positions G3 and omitting the constant 



factor 1/(2ttV0 



get from Eq. (fTS"]) : 



£ _ g- J T dt dE dc fj{E)S(t,E,c) j-jA'cv 

Jn{E)dEG x {E-E i ,5E i )^ -.8, 
S{ti^ £7, Ci c -\- Si s C(p) 



2 



(38) 

where Eq. ()35p has been rewritten using the intuitive 
shorthands cos 9 — > c, smOi — ► Sj, etc. 



IV. COMPARISON WITH THE LITERATURE 

Here, we compare our likelihood, Eq. (f3"8"|) , with cer- 
tain other likelihoods present in the recent literature and 
currently used for the analysis of SN1987A events. 

A. Jegerlehner, Neubig and Raffelt [4] 

The first likelihood is Eq. (15) of 0: 

,00 N obs 
C JNR = c - J n(E)dE -Q (3Q) 



This is an approximation of our likelihood, in that the 
background has been neglected and the time and angu- 
lar distribution are integrated (averaged) over; in other 
words, only the energy distribution is considered. How- 
ever, this expression is in direct correspondence with 
Eq. (|38|). when B — > and 86 i — > 0, and if we take as the 
definition of n(Ei) the one given in Eqs. (18,19,21) of j|. 
Furthermore, the expression of [i| agrees with Eq. 



B. Lamb and Loredo 0] 

The other likelihood that we consider is the one ad- 
vocated by Lamb and Loredo This is given by their 
Eq. (3.18) which, rewritten in our notations, reads 



C 



LL 



-J T dtdEfj(E)S(E,t) x 

UdECiiEWE^+BiEi)] 



(40) 



where we neglected dead-time, as appropriate for 
Kamiokande-II, and dropped the information about the 
angular distribution. Quoting [f|: Our derivation of the 
likelihood function reveals errors in previous attempts to 
account for the energy dependence of the efficiencies of 
the neutrino detectors; we show that these errors signifi- 
cantly corrupt previous inferences. Indeed, the likelihood 
advocated by LL has been shown to have an important 
impact for the analysis of data also by similar and inde- 
pendent analyses of SN1987A observations (l3 |. 

We would like, however, to draw the discussion on the 
correctness of the likelihood of LL. We see that the LL 
expression, Eq. (|4"0")) . coincides with our Eq. (f3"B")) only 
if we identify the function Ci with the energy response 
function of the detector: 



Ci(E) ( = } Gi(E — Ei, (Ti)r)(E), 



(41) 



where G\ is the Gaussian smearing of Eq. (|23|) . 

This is not the case for LL who instead claim (see their 
Eq. (3.21)): 



d{E) = Gi(E — Ei,ai)Q(E — E ), 



(42) 



where Eq is assumed to be the maximum energy where 
the efficiency vanishes (i.e., the minimum detectable en- 
ergy) and 6 is the step function. 

The only special case in which the LL likelihood co- 
incides with our result is when the average efficiency is 
assumed to be a step function fj(E) — Q(E — En). In 
general, this is not the case and the efficiency is a contin- 
uously growing function of the energy. The LL likelihood 
therefore incorrectly biases the analysis in favor of low en- 
ergy events. The quantitative effect of this bias on data 
analysis will be discussed further in Sect. IIVEI 

The above remarks amount to the consideration that 
the likelihood of Lamb and Loredo does not follow from 
the formal construction described in Sects. |TT] and IIIII 
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(43) 



However, it is instructive to point out more directly the 
profound principle problem of the LL likelihood. 

We begin noting that Eq. (|40|) has been derived by 
omitting constant terms from 

V LL = e - J T dtdE[f 1 (E)S(E,t)+B(E)] x 

nt c i [J dECi(E)S{E,ti) + B(Ei)] dt.dE, 

which should represents the probability that a given ex- 
perimental result is obtained. This expression is sup- 
posed to have a general validity. Then consider a simple 
limiting case: only one bin, with dimension At x AE 
and with energy above Eo; no background, B(E) = 0; 
a constant signal, S(E, t) = S; a constant efficiency, 
rj(E) = rj; a perfect energy resolution, SEi — * 0; a very 
small expected number of events, n = r/SAtAE < 1. In 
these assumptions, the most probable outcome is the case 
N cv = 0, followed by the case N cv = 1; the probability 
of other possible results is negligible. From Eq. (|43|) . we 
calculate the probabilites of the cases when no event and 
one event are observed: Vq L = 1 — n and V\ L — n/f] re- 
spectively. Their sum violates the basic principle accord- 
ing to which the sum of the probability for all possible 
results should be equal to one. 



C. Loredo and Bernstein et al. [l5| | 

We comment here on a likelihood that was not pro- 
posed for the analysis of supernova neutrinos, but that it 
is strictly connected with the previous one. 

The same position as in Eq. (|42|) was made in [3], 
where Loredo defines the quantity £i(m) — p(di\m). This 
quantity, that evidently corresponds to the quantity Ci 
discussed above, is claimed to be independent on the de- 
tection efficiency This position caused the criticism of 
Loredo [l4[ to the likelihood advocated for the analy- 
sis of trans-neptunian objects of Bernstein et al. 2004 
flil (s ee Eq. (A8) of [IH). Again, we find the position 
of [14{ unjustified, while we agree with [l5j]. In partic- 
ular, Eq. (A4) of [IH expresses the statement that the 
response function does contain the detection efficiency; it 
corresponds strictly to our Eq. (|4]). 



D. Pagliaroli et al. [13fl 

Finally, in the analyses performed by some of us [l3| , 
Eq. (|38|) was used, simplified to the case Srii — > to 
take into account the mild angular dependence of the 
v e p — > ne + reaction. This simplification does not affect 
significantly the analysis of SN1987A events. 



E. Numerical comparison of the likelihoods 

The bias implied by the likelihood of ref. [5j is nu- 
merically important already for the analysis of SN1987A 
events, as found in [5j and confirmed by [13j |. 



E 
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60 
50 



a: 40 
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i 1 1 — 

Lamb-Loredo ref.[5] 
With biasref.[13] 
Unbiased ref.[1 3] 




3 3.5 4 4.5 5 
T c [MeV] 

FIG. 2: Comparison of the joined 68% C.L. (2 dof) regions for 
the cooling radius R c and the initial temperature T c obtained 
from two analyses of SN1987A data based on the exponential 
cooling model. When Eq. (|42|l is adopted, there is a good 
agreement between the results of [f| (dotted line) and those 
of [131 ] (dashed line); the small discrepancies can be ascribed 
to the different statistical procedures (Bayesian in [5| and fre- 
quentist in [h|) and to a different numerical treatment of the 
data. The effect of switching from Eq. (|42p (dashed line) to 
Eq. (|4ip (continuous line) is much more significant [l3| . 



In order to illustrate this point better, we recall certain 
results obtained in the previous analyses. Let us begin by 
considering the conventional exponential cooling model, 
in which the v e temperature decreases exponentially with 
the time and the neutrino-radius R c is constant. As evi- 
dent from Fig. [21 the use of Eq. (gU) rather than Eq. (gT]) 
leads to important differences on the inferred values of 
the parameters. We note in particular that the (well- 
known) difference between R c and the expected size of 
the neutron star radius, R ns ~ 15 km, is amplified when 
we adopt Eq. ([42]) (i.e., when we bias the analysis). This 
outcome can be easily understood: the bias in favor of 
low energy events implies that T c (that is proportional 
to the average energy of the electron antincutrinos) will 
decrease; thus, R c has to increase to keep the number of 
events constant. 

The bias will be even more important for the analysis 
of a future galactic supernova, since the number of col- 
lected events will be much larger and the errors on the 
parameters are expected to scale as the square root of 
the number of the events. 

Indeed, the analysis of the 29 events collected by 
Kamiokande-II, 1MB and Baksan in an extension of the 
exponential cooling model leads to R c = 16^ km and 

T c = 4.6t°;6 MeV 0- A recent analysis of simu- 
lated events from a future supernova, that assumed the 
same antineutrino emission model (that includes an ini- 
tial phase of intense emission), the same central values 
as found from SN1987A (in particular, R c = 16 km and 
T c = 4.6 MeV), and a supernova located at a distance of 
20 kpc (i.e., a data set 30 times larger), leads to the con- 
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elusions that the parameters are correctly reconstructed 
when we use Eq. (|41[) . Moreover, when we combine the 
results of the simulations, we can estimate the average 
values of the parameters and of their expected errors: 
R c = 15.4 ± 0.9 km and T c = 4.6 ± 0.1 MeV [H]. 

The comparison with the values from SN1987A reveals 
that the errors are expected to decrease by about six 
times, which is similar to the improvement that we can 
ascribe to the increased number of data. We are lead 
to the conclusion that, after a future galactic supernova, 
the allowed regions in Fig. [3 should shrink by a similar 
factor in linear scale, making the effect of the bias much 
more important. 

V. SUMMARY 

We constructed the general likelihood for supernova 
data analysis, Eq. (fl5|) . and specified it to the analysis 
of SN1987A, Eq. §E§. We have compared this likelihood 



with other forms advocated in the scientific literature. 
While our likelihood is a generalization of the likelihoods 
traditionally adopted for the analysis of SN1987A events 
(or in general for the study of rare processes), it is in 
disagreement with other ones. Reasons and consequences 
of these disagreements are discussed. 
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